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Abstract 

A recently proposed reference potential approach to the inverse Schròdinger problem is further 
developed. As previously, theoretical developments are demonstrated on example of diatomic xenon 
molecule in its ground electrònic state. An exactly solvable reference potential for this quantum 
system is used, which enables to solve the related energy eigenvalue problem exactly. Moreover, 
the full energy dependence of the phase shift can also be calculated analytically, and as a particular 
result, full agreement with Levinson theorem has been achieved and explicitly demonstrated. In 
principle, this important spectral information can be reused to calculate an improved potential for 
the system, and such possibilities are discussed in the paper. Aiming at this goal, one may calculate 
an auxiliary potential with no bound states, whose spectral density for positive energies is exactly 
the same as that of the reference potential. To this end, one may solve Krein equation, which 
in the present context is simpler than using Gelfand-Levitan method. General solution of Krein 
equation can be expressed as a Neumann series. Convergence of this series of multi-dimensional 
integrals at distances not close to the origin is hard to achieve without a simple asymptotic formula 
for calculating the kernel of Krein equation. As proven in this paper, such an asymptotic formula 
exists, and its parameters can be easily ascertained. 
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I. INTRODUCTION 

Celebrated Levinson theorem 

na 

í,(0)-í,(oo)=n,7r, (1) 

originally associated with spherically symmetric potentials in three dimensions, is simple in 

form, but fundamental in content. Indeed, Eq. (1) reveals a correlation between the main 

characteristics of the continuous and discrete energy spectra, the phase shift 5i(E) (E > 0) 

and the total number (n) of bound states (E < 0). The quantity rii = n + 1/2, if the angular 

momentum l = and a so-called half-bound state exists at E = 0, and n\ = n in any other 

case. As can be shown, the solution of the Schròdinger equation for l = and E = (if 

one exists) is not normalizable, but it still gives contribution to Levinson theorem, a half 

of that compared with any "normal" bound state. Therefore, the term "half-bound state" 

introduced by Newton 3] is fitting and illustrative for this rather special solution. 

A lot of proofs of the Levinson theorem have been given by several authors using different 

methods (see, e.g., for a brief overview), and various generalizations of Eq. (1) are 

tnown. For example, Levinson theorem has been adiusted to noncentral |3f and nonlocal 
1 

5|] potentials (including bound states with E > 0), electron-atom scattering p, relativistic 
equations in one 3, two jj] and more ^ dimensions, etc. In spite of the vast literature 
on the subject, one hardly can find any examples of calculating the full energy dependence 
(from to oo) of the phase shift, which would explicitly demonstrate the Levinson theorem. 
Such a computation for a particular quantum system is a challenge even nowadays, but as 
we demonstrate below, the computational-technical dimculties can be overcome, if one is 
able to construct a reasonable reference potential for the system. Moreover, for any E > 0, 
the phase shift can be ascertained exactly, if one uses an exactly solvable reference potential. 
In the present context "exact solubility" means that the problem can be solved with the 
help of solely analytic tools, without applying methods of numerical integration. 

One-dimensional inverse theory (including Levinson theorem) can be applied to diatomic 
molècules, since a two-particle problem can be always reduced to a one-particle problem 
in an effective spherically symmetric field where the angular momentum l is replaced by 
the rotational quantum number J. On this basis, the reference potential method developed 
in this paper is illustrated on example of diatomic xenon molecule (Xea) in its ground 
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electrònic state. To avoid excessive technical difnculties and concentrate on the essence of 
the approach, only the rotationless case (J = 0) is analyzed. 

Complete knowledge of the phase shift is one of the three obligatory preconditions for 
the unique solution of the inverse Schròdinger problem for a quantum system. In addition, 
one has to ascertain full energy spectrum of the bound states E m < (m = 1, ...,n) and n 
associated real parameters C m that uniquely fix their normalization lOj . In principle, all 
these data can be obtained experiment ally, but unfortunately, this is almost unachievable 
for a real system. In this connection, one may ask: is it still possible to get some practical 
benefit from the rigorous methods of the inverse quantum theory? One option, as already 
mentioned, is to construct a reference potential based on the known data on the system's 
properties. Since the reference potential is fixed, all its spectral parameters are determined a 
priori. Naturally, there is no need to regain a potential which is already known by defmition, 
but the idea is to use the calculated reference parameters as zeroth approximations to the 
real characteristics. Using these approximate spectral characteristics, it is possible, at least 
in principle, to build up a more rea 



construct a Bargmann potential 



11 



istic potential for the system. For example, one can 
, Q whose bound states coincide with the actually 
observed discrete energy levels, or adjust the potential more adequately to the available 
scattering data. These possibilities will be discussed in this paper, which is organized as 
follows. 

In Section II we first specify the analytic form of the reference potential, and then describe 
the details of calculating the phase shift for the full range of scattering states. Having 
ascertained the phase shift, one can calculate another important spectral characteristic, Jost 
function, and then fix spectral density (the terms to be specified below) for the reference 
potential. In Section III we demonstrate how these spectral characteristics can be reused 
to calculate a special auxiliary potential with no bound states, applying Krein method 
of solving the inverse scattering problem. Starting from this auxiliary potential, one can 
calculate an improved (Bargmann) potential for the system, as already mentioned. Several 
possibilities of this kind are discussed in Section IV. Finally, the main results of the work 
are summarized in Section V. 
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II. FULL ENERGY DEPENDENCE OF THE PHASE SHIFT 



According to the scheme outlined in the previous section, our first goal is to construct 
an exactly solvable (in the above specified sense) reference potential for the system. A 
good choice for this purpose, as explained in detail elsewhere is a multi-component 

potential composed of several smoothly joined Morse type pieces 

V(r) = V k + D k [exp(-a fe (r - r k )) -lf,r6 (0, oo), (2) 

where V k , D k , a k and r k are some constants (not definitely positive), and the subscript 
k corresponds to different components smoothly joined at some suitably chosen boundary 
points X k+ i. A five-component reference potential for Xe2 has been constructed |14| . which 
is in good agreement with all known spectroscopic data. In this paper the emphasis is 
put on principies, not on specific details. Therefore, a more simple (but still reasonable) 
reference potential of three Morse components (k = 0,1,2) is used. The most internal 
of them (k = 0) represents a so-called pseudo-Morse potential. It means that the tiny 
potential well corresponding to this component (if taken separately) is just of the limit 
depth to completely lose the discrete energy spectrum. Consequently, D = h 2 al/(8m) (m 
being the reduced mass of the pair of àtoms), so that only three independent parameters 
remain. The central component (k = 1) is an ordinary Morse potential, while the most 
external one (k = 2) is a "reversed" Morse potential with the parameter D2 being negative. 
By introducing a "reversed" component one artificially creates a potential barrier in the 
long-distance range. This might seem unphysical and unjustified, but the point is that 
the height of this artificial hump approaches zero as the parameter r 2 approaches infinity. 
Therefore, taking a sufficiently large r 2 , the hump becomes almost insignificant, while the 
analytic treatment remains simple and flexible. Almost equivalently, one may fix V(r) = 0, 
if r > r 2 . The parameters of the reference potential are given elsewhere |l5j |. and the curve 
itself is shown in Fig. 1. The described construction based on Morse type functions is not 
occasional. On the contrary, this is one of the very few options to easily get two linearly 
independent analytic solutions of the Schrodinger equation for the whole physical domain. 

As already mentioned, the phase shift for the specified reference potential for any 
E e (0, 00) can be ascertained analytically. To demonstrate this, let us briefly recali the so- 
lution scheme for Morse- type potentials (see for more details). First one introduces 

dimensionless variables y k = 2a k exp(—a k (r — r k )), where the subscripts k = 0,1,2 have 



the same meaning as specified above, and a k = \/2mD k / (Ha k )- The relevant Schròdinger 
equations then read 



d 2 V(y k ) 1 dt>(y k ) 



dy 2 k ' Vk dy k 



j4 ± (°± _ 1 

Vi \Vk 4 



o. 



(3) 



where the plus sign in square brackets corresponds to the subscripts k = and k = 1, and 
minus - to k = 2. The quantities n\ are defined as follows: n\ = (a\/D k ) ■ (Vk + D k — E). 
Note that D 2 < and ao = 1/2. Next, Eq. (3) is converted to confluent hypergeo- 
metric equation, whose two linearly independent solutions can be always ascertained an- 
alytically. For example, general solution for the internal region r < X\ can be writ- 



ten as #„(0 = A^J^r) + N 2 ^'(r), where ^'(r) = A Q (y ) cos [B Q (y ) - ip Q - a ol· i r], 
A (y ) sin [B (y ) — <p — a /i r], the functions A (y ) and B (y ) are defined as 



,(2) 



tt«(r) 
follows: 



^od/o)e iBo(w) = exp(-yo/2)iF 1 (z// ; 2ifi + 1; y ) 



(4) 



and i Fi fa; c; x) denotes Kummer confluent hypergeometric function. As demonstrated else- 
where special solution (r) is almost insignificant if a classically forbidden region 
exists, i.e., if E < V(0) (note that V(0) is very large, but still finite in the framework of 
our approach). Correspondingly, the correct physical solution (apart from normalization) 
practically reduces to the particular solution \l/[ ) 1 ' ) (r).Therefore, the phase parameter (NB! 

r tv?.?' i i^i 

for the region E < V(0) can be 

_ i i //ió) j 

calculated very easily 



A*o ( "o^o + 1 - ln 2 - - ln(l + 4# 



+ 



oo 

\I( 



1 \ . . .dt 
( othí — — 1 e sin(2/i í) — . 



(5) 



where the integral can be conveniently evaluated 

T 

r í 1 \ d* 

7 

o ' o 







oo T 

y (cothí - e- í sin(2^)y = ^ e"' sm(ir^)f(t)dt, T = ^ 



(6) 



/(*) 



coth í — T 
t 



-T 



coth(t + T) - -L 



t + T 



^ _ 2T coth(t + 2T) 



i 



í+2T 



Another equivalent formula for this quantity reads 



t + 2T 



n=l 



m— lo2n 



7 n 



2 2n 5 r , 



(2n)(2n-l)(l + 4^g 



2N2n- 



fc=0 V 7 



(7) 



Here B n denotes the n-th order Bernoulli number. Since two different definitions of the 

Bemoulli a Umber8 a r e w.de.v ased Q (and —es co^ed), we have to 8 pec lfy th- 

the "old-style" definition is meant here, i.e. B% — 1/6, _B 2 = 1/30, B 3 = 1/42, etc. 

The simple scheme just described cannot be used in the region E > V(0) where special 

solutions (r) and \Pq (r) become equally important. Nevertheless, the analytic treatment 

can be preserved. Indeed, phase shift is related to regular solutions of the Schròdinger 

equation, which means that the physically correct linear combination of v&q (r) and \&q (r) 

should vanish as r — > 0. From this one immediately concludes that the correct wave function 

\l/(r) ~ sin{a /io^ + B [yo(0)] — B [yo (?")]}■ According to Eq. (4) it means that one first 

has to ascertain the argument of the complex function xF\ [i/J-o] 2i/^o + 1; 2/o( r )] at r = 0, to 

calculate the phase shift for the region E > V(0). This is a really complicated task, if one 

tries to use primary definition 

„ . , a z a(a + 1) z 2 

1 F l(o;c;2; )^l + - ïï + -i_ T 2- + ...., (8) 

because z = yo(0) is a very large quantity, and as a result, the series converges very-very 
slowly. It means that extremely high precision of computations is needed to get the correct 
result. The convergence can be essentially improved due to the relationship between the 
arguments, c = 2a + 1, specific to the pseudo-Morse approximation [14], namely 

( 9 ) 

where (a) n = T(a + n)/T(a) is the Pochhammer symbol. There exists even third, and in 
the present context, most appropriate alternative to calculate the needed quantity Bq(z). 
Namely, iFi(a;c;z) can be expressed in terms of Buchholz polynomials [17| Pn(b,z), and 
as has been shown these polynomials in turn can be given as sums of products of 
polynomials in b (related to E) and z (related to r), separately. In the present case this 
leads to the formula 

^=È fc (%H <fc'} , (10) 

oo z n 

where o-^iC a; z) = ^2 , and the Buchholz polynomials read 

n=0 n\(a) n 

s=0 \ ' 



with / (6) = 1, f s (b) = -(b/2 - 1)E(V) - ^^ frib), and g (z) = 1, (/ ro (z) = 

r=0 s — r 

j z I 2 ^] 4* +1 .B fc+1 

_ T S ("fc ) ,1-, — 9m-2k-i{z) (s = 1,2,..., and m = 1,2,...). Again, as above, fi n 
denotes the n-th order "old-style" Bernoulli number. 

Thus we have described all details of calculating the phase parameter within the pseudo- 
Morse approximation, which is used for r G [0, . Having fixed this parameter, one can eas- 
ily calculate the real phase shift, using the general asymptotic formula \l/(r) ~ sin [kr + S(k)] 



(k = \J 2mE / h) , and the continuity conditions for the wave function and its derivative at 
the boundary points Xk+i- The whole procedure is based on exact analytic formulae, and 
there is no need for numerical solution of the well-known phase equation [19] 

S'( r , k) = _^mV(r) ^ + ^ ^ = ^ 

nk 

which would be a much more difficult task, since the phase shift can only be obtained as 
the limit S(k) = lim 5(r,k). We still have to mention that it is not so easy to make a 

r— >oo 

reasonable use of Eqs. (9) or (10), in order to ascertain the quantity B [yo(0)] . Although 
convergence of these complex series is much better compared with Eq. (8), one definitely has 
to overcome the barrier of Standard double (or fourfold) precision of high level programming 
languages, because much higher accuracy is actually needed. Fortunately, this is not any 
serious problem nowadays. Fast and powerful computational tools are freely available, such 
as high-precision environment UBasic and an arbitrary precision library MAPM for 
C/C++ [21]. These tools have been used to ascertain the energy dependence shown in Fig. 
2. One can see a really good agreement with the Levinson theorem, but to achieve this, 
the phase shift has been calculated up to very high energies (note that the energy scale in 
Fig. 2 is logarithmic, and it involves 20 orders of magnitude!). A good point is that there is 
no need for cumbersome high-precision calculations at very high energies, because this part 
of the curve can be accurately ascertained from general considerations. Namely, from Eq. 

i oo fe2 

(12), one gets S(k) = lim S(r, k) = — — j V(r) {1 - cos [2kr + 25{r, k)}} dr (C = ). 

r— >oo 2Ck o 2m 

Repeatedly integrating by parts and assuming k — > oo, the following asymptotic formula 



can be obtained 



sí - k) = T + ¥ + ¥ + -· k ^ OD ' (13) 



f V(r)dr 



where ai = — — and other coeficients 03, a 5 ,... are also directly related to the 
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potential and its derivatives (see for more details). Using Eq. (13) one can extend 
calculations to arbitrarily high energies, and thus, as a particular result, explicitly verify 
the Levinson theorem. Since the validity of Levinson theorem is beyond any doubts, Fig. 2 
actually demonstrates the success of the described approach to calculating the phase shift. 

It might be of interest to compare the exact phase shift with the relevant quantity obtained 
from a semi-classical (WKB) calculation. In the framework of the WKB approximation, the 



phase shift reads 



22 



S(E) 



ro 



E - V{r) 1 



C 



7T 



dr + — — kr , 



(14) 



where k 



p! and r is the solution of the equation 



E - V(r) 1 



——. In Fig. 3 one can 

4H 



c — u " C 

see a part of the exact phase curve in comparison with the WKB phase shift, calculated 
according to Eq. (14). As expected, semi-classical model becomes more and more reliable 
with increasing energy, but it fails to adequately describe the energy range E < V(0) where 
the exact phase can be ascertained very easily. 



III. PREPARATORY STEPS TO APPLY KREIN METHOD 

How could one get any practical benefit from full knowledge of the phase shift and 
other spectral characteristics related to the reference potential? Trying to answer this 
question, let us introduce another important spectral characteristic, Jost function F(k) = 
|-F(fc)| exp [— i5(k)], which fixes relationship between regular and physical solutions of the 



Schròdinger equation: \I/(r, k) 



k 



F(k) 

F(E)\ = ]J(1-E n /E) exp 



N 



ip(r, k). The modulus of Jost function reads 
8{E')dE 



n=0 



1p 

lï 



E'-E 



E e (0,oo). 



(15) 



Here E n are the discrete energy eigenvalues and symbol P denotes the principal value of the 
integral. Next, we can fix spectral density 



dp(E) 
dE 



ir-WE\F{E)\~ Z ,E > 0, 
J2C n 5(E-E n ), E<0. 



(16) 



which is closely related to the real solution schemes of the inverse problem. An interesting 
option, among others, is to calculate an auxiliary potential Vb(r) with no bound states, 
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whose spectral density for E > is exactly the same as that of the reference potential 
(\F(E) \ remains the same). We first show how this problem could be solved and then, in the 



next section, discuss of the usefulness o 
Vo(r), one can apply Krein method 23, 



the approach for practical purposes. To calculate 



24 1, based on the following integral equation: 



2r 



r 2r (r') + H{r') + J T 2r {s)H{s - r')ds = 



whose kernel reads 



H(r) 



7Ï 



-1 



— 1 cos(kr)dk. 



The desired potential with no bound states then becomes [ï(| 



(17) 



(18) 



V (r)=4C^[G(x 



íi 



dG(x) 
dx 



, x = 2r, G(x) = r 2r (2r). 



(19) 



Since Eq. (16) is a Fredholm integral equation of second kind, its solution can be given 
Neumann series 



G(x) = —H(x) + / H(x — x\)H{x])dx\ —II H{x — x 2 )H(x2 — X\)H (x\)dx2dxi + . 




o oo 

XX X 



+(-l) n J j ... J H(x - x n )H(x n - x n -x)...H(x2 - x 1 )H(x 1 )dx n dx n - 1 ...dx 1 + ... (20) 



This general formula can be easily proved, if one uses a suitable quadrature rule to discreticize 
Eq. (17), then finds the last element of the solution vector (the only one which is actually 
needed) of the corresponding matrix equation, and finally goes to the limit h — *■ 0, where h 
is the step of integration. For relativel y sm all distances the discrete analogue of Eq. (17) 



>r relativel y si 
imination 26j 



can be solved by Gaussian elimination [26j, but with increasing x the number of equations 
rapidly becomes too large for this kind of approach. In this situation, one is guided to 
multidimensional Monte Cario techniques using, for example, Metròpolis algorithm j^. 
This interesting possibility, however, is a subject for a forthcoming study and will be not 
further discussed here. 

Calculation of Jost function and the related function q(k) = ^ — 1 is described 

nn imi 

elsewhere |l5l.l26j|. As can be seen in Fig. 4, g(k) = —1 (with high accuracy) in a wide range 



below (and not too close) to the characteristic wave number k = ^2mV(0)/h, and accurate 



analytic approximations can be obtained for the remaining part of the curve. Again, as for 
the phase shift, a simple asymptotic formula 

^ ) = F + F + F' k ^ K (21) 

can be derived, where the leading coefficient b\ = — ^77^ an d the value k a = 75000 A _1 
has been fixed for the model system (Xe2) studied. Since H(r) = H(—r) is just the Fourier 
cosine transform of the characteristic function g(k), calculation of this kernel in the short- 
distance range can be performed quite easily and accurately. As shown in Fig. 5, H(r) is 
a rapidly oscillating function with decaying amplitude. At longer distances, however, when 
the amplitude of H(r) becomes very small, even minor inaccuracies in calculating g(k) result 
in large relative errors of its Fourier transform. Fortunately, a correct asymptotic formula 
for H(r) can be fixed, namely 

H(r) a exp(— br) cos(fcr + a), r —>■ 00, (22) 

where a,b,k and a are some constants to be specified. To prové Eq. (22), let us assume 

bi 

that k — > 00, so that according to Eq. (21) one can approximate g(k) = — . We can 

k 2 

00 

perform the inverse Fourier transform g(k) = 2 J H(r) cos(kr)dr, and divide the domain into 



two intervals, (0, r a ] and [r a , 00), where r a is an arbitrary (but sufficiently large) quantity, 
such that sin(/cr a ) = and cos(/cr a ) = 1 (these conditions can be always set, because 
k is arbitrary). Then, taking account of the general formula H'(0) = 



2 4C 

integrating by parts in the first interval, and using Eq. (22) for the second interval, one 

00 

comes to a relation H'(r a ) + ak 2 j exp(-br) cos(kr + a) cos(kr)dr = 0. It can be explicitly 
checked that this relation holds as k —>■ 00, which completes the proof. 

Since |-F(0)| is an immensely large quantity Q], g(0) = 2 J H(r)dr = —1 (with very 



high precision). Let x be an arbitrary distance, such that Eq. (22) can be used to calculate 
H(x). Then 

X 00 

J H{r)dr + 1/2 = - J H(r)dr = ^^"^ cos{kx + a + f3), (23) 

X 

b k 

where cos a = _ : and sin a = —====. Thereafter, accurately evaluating the integral 

Vb 2 + k 2 Vb 2 + k 2 

X 

j H(r)dr, the desired parameters can be easily fixed by direct least-squares fit of Eq. (23). 
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To be more specific, Simpson rule with integration step h = 10 7 A was used for x = 0.038 
A, and as a result, the following set of parameters has been ascertained: a = 2.7125 1/À, 
b = 70.04656 1/À, a = -3.67004, and most importantly, k = 19165.07 1/À ta 0.99866 
ho- Fig. 6 demonstrates how nicely Eq. (22) already works at r > 0.03 A. Thus, for this 
region there is no need for troublesome Fourier transform to calculate the kernel of the Krein 
equation, which makes Eq. (20) much more useful for practical purposes. 



IV. IMPROVING THE INITIAL POTENTIAL 



Having ascertained all spectral characteristics of the reference potential and described a 
possible solution scheme of the Krein equation, one cannot ignore the principle question: 
what is the practical point of all this? Can these data be used to calculate an improved 
potential for the system? First, we have to point out that such an improved (in some sense) 
potential V\{r) can be ascertained without auxiliary calculations described in Section III. 
Most easily this can be demonstrated in the frame of Gelfand-Levitan method (2^, which 
is based on the integral equation 

r 

K(r, r') + G(r, r') + J K(r, s)G(s, r')ds = 0, (24) 


whose kernel can be directly linked with the reference potential V(r) 

00 

G(r,r>) = j [dp {1 \E)-dp(°\E)] tp(r, E)ip(r',E). (25) 

—00 

Here 

dp^(E)-dpW(E) = í *- l ^E [\FW(E)\- 2 - \F(E)\- 2 ] , E > 0, 

dE \Y.c£ ) 8{e-eV)-y<c% ) 8{e-e ( $), e<o, i2b) 

\ k m 

the regular solution ip(r,E) as well as the quantities p^(E), \F( \E)\ , Em and Cm\ 
(m = 1,2, ...,n) correspond to the reference potential V(r), while dp^(E), \F^(E) \ , E^ 
and Cfr (k — 1 ; 2, ...,/) are related to Vi(r). Suppose we are looking for a potential Ví(r) 
whose spectral density at E > is exactly the same as that for the reference potential 
V(r), i.e., \F^(E)\ = \F(E)\, but the discrete energy spectrum and the related norming 
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constants are different. In this case the kernel is separable: 

G(r, r>) = £ CÏ% (r, E?) <p (r', E?) - £ £7£V (r, üff) ? (r', , (27) 

fe m 

and Eq. (24) can be easily solved. Thus, in principle, one can directly calculate the difference 
of the two potentials in question (see, e.g., [ï(| for more details) 

Vx{r)- V{r) = 2C-^-K(r,r) = -2C {ln [det S (n+/) (r)] }" , (28) 

r 

where the (n + l) x (n + l) matrix S^ n+l '(r) = I + f R^ n+l \s)ds, I is the corresponding unit 

o 

matrix, and the elements of another matrix R^ n+l '(s) read 

Rf k +l \s) = C jV > (s, E 3 ) <p (s, E k ) , (j, k = l,2,...,n + /). (29) 

The quantities Cj and Ej are defined as follows: Cj = Cf\ and Ej = Ej, if 1 < j < l; 

Cj = — Cj_j, and £/j = if l < j < n + l.ïï the parameters iíj and Cj are taken equal 

to the actually observed vàlues of these quantities, the potential V\{r) can be considered more 

realistic than the reference potential V(r). 

According to Eqs. (28) and (29), one only has to accurately evaluate the regular solutions 

íç (s, Ej) (j = 1, 2, ...,n + l), in order to ascertain the improved potential Vi(r), which seems 

to essentially reduce the computation time. This imagination, however, is deceptive, because 

full knowledge (from to oo) of n + 1 regular solutions is in fact equivalent to calculating 

n + l auxiliary potentials for the system. Indeed, let us recali that once the ground state 

wave function \&o( r ) (or any other bound state wave function) is known, the potential U(r) 

is immediately known as well (up to a constant) j^J. This general rule directly results from 

the Schròdinger equation, if one takes the ground state energy as the zero point of the energy 

scale: U(r) = C \E'g(r)/\í'o( r )· I n the present case we have to take into consideration that 

the solutions (p(r,Ej^) are not the real eigenfunctions, if Ej ^ Nevertheless, since 

the norming constants C^ are expected to be known, the regular eigenfunctions <fi(r, E^) 

related to the potential VAr) can be easily ascertained using the well-known relations of the 

quantum inverse theory pjj] . For example, if there was only one bound state with energy 

r r i-i 

E^and norming constant C-[ , then (pi(r,Ef') = (p(r,Ef*) 1 + C± J <p 2 (s, E^)ds 

L o 

Thus, seemingly simple Eq. (28) "hiddenly" contains a procedure of removing all previous 
bound states (equivalent to calculating regular eigenfunctions tp(r,Ej), j = 1,2, ...,n ), 
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and introducing new bound states (equivalent to calculating regular solutions (p(r, E^), 
j = 1,2,...,/ ) related to the improved potential V\(r). 

The explícit procedure equivalent to Eq. (28) can be performed as follows 30. 
Starting from the potential V n (r) = V(r), one first removes its zeroth level Eq and calculates 
a new potential, isospectral with the initial one (apart from Eq) 

2 



K-i(r) 



V n {r) 



2C 



+ 



ds 



(30) 



Here W Q n \r) is the eigenfunction of the zeroth level, not necessarily normalized (note that 
the norming constant is absent here). The new set of wave functions for V n -i(r) can also 
be expressed in terms of the previous wave functions for V n (r) íÜj]. The procedure can be 
continued, i.e., one can remove, one by one, all bound states along with calculating new 
auxiliary potentials at each step, until he finally comes to the potential Vo(r) with no bound 



states. Thereafter, one can introduce new bound states £7- with norming constants C-^in 



exactly the same way as described above, i.e., 

Vl(t) - V (r) = -2C {ln [det S®(r)] }" 

r 

where the / x / matrix (r) = I + f R^(s)ds, and 



(31) 



*8« 



C®<p{8,EW)<p{8,EP), (j,k = l,2,...,l). 



(32) 



Alternatively, one can solve Krein equation, as described in Section III, and then use 
Eq. (19) to get the auxiliary potential Vo(r) with no bound states. As mentioned, all three 
approaches to calculating the improved potential V\{r) are equivalent, in presumption that 
the modulus of Jost function remains unchanged, which means that 1 101 




and 




■F (k), 



(33) 



(34) 



where Fi(k), F(k) and F Q (k) are the Jost functions for the potentials Vi(r), V(r) and Vb(r) 



respectively, k) 



2mE) '/n and k) 



2mEf ) /h 
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An alternative to applying Eq. (31) is to introduce new bound states step-by-step, start- 
ing from the auxiliary potential Vo(r). This well-known procedure is based on the formula 



10| (see also |26|) 



V k (r)=V k - 1 {r)-2C{hi 



r 



(k = l,2,...,l) 



(35) 



where the regular solution ^_ 1 (s, Ejp) is not the real eigenfunction, since it is related to 
the potential V k -i(r). Naturally, the norming constants Cf ] in Eqs. (29), (32) and (35) 
are different (although the same notation is used), because they are related to different 
potentials. 

Currently, one can only guess which of the described approaches is computationally most 
effective. However, one cannot argue that solution of Krein equation is an interesting prob- 
lem on its own, and the scheme outlined in Section III is realizable. At least, it seems 
reasonable to first calculate the auxiliary potential Vq(t) with no bound states, and then try 
to improve the initial potential. Introducing "real" bound states is not the only option in 
this context. For example, one can construct another potential Vo(r) with no bound states, 
a Bargmann potential whose Jost function reads 

' k + iüj 



F (k) 



n 

i=i 



k + ibj 



F (k), 



(36) 



where the quantities aj ^ ò,- are real and positive. For simplicity, let us assume that m = 1 



y 



(ai = a and b\ = b). Then 

V (r) - V (r) 



-2C l ln 



b 2 -a? J 



(37) 



where the symbol W denotes Wronskian determinant, and fi(ia,r) is the Jost solution of 
the Schròdinger equation (fi(r,ia) — > exp(— ar) as r — > oo). Both f\(r,ia) and the regular 
solution <fi(r,ib) are related to the potential Vo(r). Thus, suitably choosing the parameters 
cij and bj, one can, in principle, take a more adequate account of the scattering data or of 
the known behavior of the potential near the origin. Suppose, for example, that the quantity 
Vq(0) is known. Then, using the formula 



d 
dr 



W [/i(r,m),(ii(r,ift)] 
b 2 - a? 



-fi{r,ià),<pi{r,ib), 



(38) 
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which is vàlid for any two solutions of the Schròdinger equation for the same potential, 



7(0) _ - 

taking account that -Fb(^) = 1 + - + 0(k ) 15j . and consequently, -fo(^) = 1 + 

26 k 2 

7(0) 



2CP 



2Cfc 2 

0(A;~ 4 ), one comes to a more transparent form of Eq. (37): 



7o(r) - 7 (r) _ d 
7 (0) - 7 (0) rfr 



/i(r,zo),^i(r,i6) 



.Wl/ifaiaJ^r.ifc)]. • (39) 
One can easily check that the right side of Eq. (39) indeed approaches unity as r — > 0, 
because <fi(0, ib) = 0, but /i(0,za) ^ (otherwise E a = —Ca 2 would be a bound state). 

Using the modified auxiliary potential 7o(r), one can construct another improved po- 
tential for the system, introducing new bound states all at once, according to Eq. (31) or 
step-by-step, according to Eq. (35). 

V. CONCLUSION 

Starting from a known reference potential, one can calculate important spectral charac- 
teristics of the quantum system, which are almost unattainable in a real experiment. Then 
it is possible, in principle, to construct another, in some sense more realistic potential for the 
system. This kind of approach has been applied to diatomic xenon molecule in its ground 
electrònic state. The phase shift for the full range of scattering states has been calculated, 
and excel·lent agreement with Levinson theorem explicitly demonstrated. All this impor- 
tant spectral information, along with full discrete energy spectrum, can be stored into Jost 
function, and then reused for constructing a new (possibly improved) potential. 

Choice of a proper reference potential, applicable in the whole physical domain, is a 
difficult and rather speculative task in any particular case. A variety of semiclassical schemes 
have been devised to deduce the potential directly from the available experimental data, and 



these concepts are constantly improved and developed (see, e.g., J32j,|33|]). Such methods, 
although inaccurate from the rigorous quantum mechanical point of view, have proven useful 
for practical purposes. Therefore, they might help to construct a reference potential, which 
can be used as a reasonable approximation to the real potential. Only if the reference 
potential is realistic enough, there is a chance to construct a Bargmann potential that would 
be even more realistic. On the other hand, as demonstrated in this paper, the còmput ational- 
technical difficulties can be overcome much more easily, if one chooses an exactly solvable 
reference potential (in the sense specified in Section I). To this end, a multi-component 
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potential of smoothly joined Morse-type pieces is especially suitable. Full energy dependence 
of the phase shift for such a potential can be ascertained analytically, using a method 
(described in Section II) , which is largely based on a so-called pseudo-Morse approximation 
for the small distances (and high energies) region of the potential. In this paper, only a 



single pseudo- 
elsewhere 



VIorse component has been used for the range r £ [0, Xi\. As demonstrated 



141 ] . adding more pseudo-Morse components does not bring along any serious 



computational problems, so one can include just as many components of this type as he 
considers reasonable. 

An important constituent of the overall concept is Krein method for calculating an aux- 
iliary potential Vq(v) with no bound states. For this purpose one uses the Jost function, 
which has been ascertained for the reference potential. General solution of Krein equation, 
Eq. (20), can be given as a series of multi-dimensional integrals (Neumann series). Thus, 
Eq. (20) represents a formulation of a serious (but very interesting) computational-technical 
problem, not solved yet. Fortunately, a simple asymptotic formula, Eq. (22), for the kernel 
of Krein equation can be used to accurately evaluate the multi-dimensional integrands in 
Eq. (20). As demonstrated in Section III, Eq. (22) can already be used at rather small 
distances, r > 0.03 A, instead of performing Fourier cosine transform of the characteristic 
function g(k) there. 

We discussed several possibilities of improving the initial reference potential in the rigor- 
ous framework of the quantum-mechanical inverse theory. In particular, explicit calculation 
of the auxiliary potential Vo(r) is not obligatory for this purpose, if one prefers to calculate 
regular eigenfunctions of all bound states for the whole physical domain, and additionally 
assumes that |F^ 1 ^(_E')| = \F(E)\ = \F (E)\ (see Eqs. (33), (34) and the paragraph between 
Eqs. (26) and (27)). In this case, one can try to use Eqs. (28) and (29) to ascertain an 
improved potential V\(r) whose discrete energy spectrum would exactly coincide with the 
"real" one. This method is extremely difficult to implement, but exactly the same can be 
said about other possible approaches described and discussed in Section IV. None of these 
schemes has been realized yet, and there is no real ground to give preference to any of 
them. In this paper the emphasis has been put on promoting Krein method to calculat- 
ing the auxiliary potential with no bound states, which is an interesting theoretical and 
computational-technical problem on its own. Ràpid development of parallel computing and 
Grid technology opens a new powerful channel for the solution of such complex problems, 
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which suggest optimism for further research in this direction. 
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Figure captions 



Fig. 1. Three-component model potential for the system (Xe 2 in ground electrònic state) in- 
vestigated. All components have the well-known analytic form of the Morse potential, 
but the ordinary Morse approximation is used only in the central range r G [Xx, X 2 ] 
(see the explanations in Section II). The parameters of the components as well as the 
calculated discrete energy levels (24 in total) are given in Ref. 15. 

Fig. 2. Explicit demonstration of the Levinson theorem (5(0) — 5(oo) = nn) for the model 
system studied. As needed, 5(0) = 24ir, since the system has 24 bound states. At 
E = 3.146294 meV, the phase shift passes a zero, and then remains negative, very 
slowly approaching the limit (S(oo) = 0) as E — > oo. The left side-inset shows the 
nearly linear energy dependence as E — > 0, in full agreement with general-theoretical 
concepts. In the right-side inset one can see that the phase curve has an inflection 
point near E = V(0). 

Fig. 3. A comparison between the exact phase shift (sòlid line) and WKB phase shift (open 
circles), calculated according to Eq. (14). The difference between the two phase curves 
is shown in the inset. 

Fig. 4. Visualization of the characteristic function g(k) = \E(k)\~ 2 — 1, which determines the 
kernel of the Krein equation. As can be seen, g(k) = —1 (with very high precision) if 
k ~ = ^/2mV(0)/h, then starts to monotonously decrease, and g(k) — > as k — > oo. 
In spite of the seeming simplicity of the energy dependence, the whole curve has to 
be calculated very accurately up to high energies, in order to accurately ascertain the 
kernel H(r). The high-energy part of the curve (k > k a ) has been calculated according 
to Eq. (21). The asymptotic wave number k (see Eq. (22) ) is also shown. 

Fig. 5. Demonstration of the behavior of Krein ií-function (kernel of Krein equation) given 

by Eq. (18). The lower graphs start where the upper ones end. Note that the period 

2tt 

of oscillations rapidly achieves the characteristic value -p. 

K 

Fig. 6. Demonstration of the validity of the asymptotic formula for calculating the kernel 
of Krein equation. Open circles in the figure represent the result of Fourier cosine 
transform according to Eq. (18), while the sòlid line corresponds to Eq. (22). One can 
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see excel·lent agreement between the two principally different approaches to calculating 
the same quant ity. 
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